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Anisotropic flows (ui, ^ 2 and ^4) of light nuclear clusters are studied by a nucleonic 
transport model in intermediate energy heavy ion collisions. The number-of-nucleon scal- 
ings of the directed flow (yi) and elliptic flow (v 2 ) are demonstrated for light nuclear 
clusters. Moreover, the ratios of v±jv\ of nuclear clusters show a constant value of 1/2 
regardless of the transverse momentum. The above phenomena can be understood by the 
coalescence mechanism in nucleonic level and are worthy to be explored in experiments. 

1. Introduction 

Anisotropic flow can reflect the initial state of the heavy ion reaction and evolution 
dynamics. In addition, the information of the nuclear equation of state (EOS) could be 
also learned from the flow studies [[HElEllllElElIZllHlini HDl dl . Many studies of 
the 1-th and 2-nd anisotropic flows, namely the directed flow and elliptic flow revealed 
much rich physics in heavy ion collision (HIC) dynamics. Recently, it was demonstrated 
that the number of constituent-quark (NCQ) scaling exhibits from the measurement of 
transverse momentum dependence of the elliptic flow for different mesons and baryons 
in ultra-relativistic Au + Au collision experiments at the Relativistic Heavy Ion Collider 
(RHIC) in Brookhaven National Laboratory [021? it indicates that the partonic degree of 
the freedom plays a dominant role in the formation of the dense matter during early stage 
of collisions. The coalescence or recombination models of the constituent quarks have been 
proposed to interpret the NCQ-scaling of hadrons at RHIC [ O HI HSl HEl H3 HH| • On 
the other hand, the coalescence mechanism has been also used to explain the formation of 
light fragments and their spectra of kinetic energy or momentum in intermediate energy 
heavy ion collisions some times ago [O (201 EH E21 E>j In addition, mass dependence 
of the directed flow was experimentally investigated in a few studies [ [21 However, 
systematic theoretical studies on the anisotropic flow of different nuclear fragments in 
intermediate energy domain in terms of the coalescence mechanism is still rare. 

In this paper, we will investigate the flow scaling of light nuclear fragments in a frame- 
work of a nucleonic transport model, namely isospin-dependent quantum molecular dy- 
namics (IDQMD) at intermediate energies. A naive coalescence mechanism in nucleonic 
level was applied to analyze the fragment flows. Dependences of the number-of-nucleon 
(NN) of the anisotropic flows, v\ and v 2) are surveyed, and the ratio of V4 to v\ is studied. 
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2. Definition of anisotropic flow and Model Introduction 

2.1. Definition of anisotropic flow 

Anisotropic flow is defined as the different n-th harmonic coefficient v n of an azimuthal 
Fourier expansion of the particle invariant distribution [ Ej 

dN 00 

— - oc 1 + 2 J2 v n cos(n4>), (1) 

# n=l 

where is the azimuthal angle between the transverse momentum of the particle and the 
reaction plane. Note that in the coordinate system the z-axis along the beam axis, and 
the impact parameter axis is labelled as x-axis. 

The first harmonic coefficient Vi represents the directed flow, 

v x = (cos® = <^>, (2) 

Pt 



where p t = ^p\ + p 2 y is transverse momentum. v 2 represents the elliptic flow which 
characterizes the eccentricity of the particle distribution in momentum space, 

v 2 = {coaW)) = (3) 

Pt 

and represents the 4-th momentum anisotropy, 

/ Pt ~ Qp 2 X P 2 y + Py\ 

l " = \ — Vi — /■ (4) 

In this work, we will study the above flow components. 

2.2. Nucleonic transport model: QMD Model 

The Quantum Molecular Dynamics (QMD) approach is an n-body theory to describe 
heavy ion reactions from intermediate energy to 2 A GeV. It includes several important 
parts: the initialization of the target and the projectile nucleons, the propagation of 
nucleons in the effective potential, the collisions between the nucleons, the Pauli blocking 
effect and the numerical tests. A general review about QMD model can be found in [ 26j. 
The IDQMD model is based on QMD model affiliating the isospin factors, which includes 
the mean field, two-body nucleon-nucleon (NN) collisions and Pauli blocking [ E3 EHl ESj • 

In the QMD model each nucleon is represented by a Gaussian wave packet with a width 
\/T (here L = 2.16 fm 2 ) centered around the mean position fl(t) and the mean momentum 

( ^ = ( 27r L) 3/4 exp[ iZ — ]exp[ n — ] * ® 

The nucleons interact via nuclear mean field and nucleon-nucleon collision. The nuclear 
mean field can be parameterized by 

U(p, Tz ) = <*(-£-) + (3{^r + r z )V c + C sy J Pn - Pp) r z + U Yuk (6) 

Po Po 2 p 
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with p the normal nuclear matter density (0.16 fm~ 3 is used here), p, p n and p p are the 
total, neutron and proton densities, respectively. r z is zth component of the isospin degree 
of freedom, which equals 1 or -1 for neutrons or protons, respectively. The coefficients a, 
(3 and 7 are parameters for nuclear equation of state. Two set parameters are used: a 
= -124 MeV, (3 = 70.5 MeV and 7 = 2.0 which corresponds to the so-called hard EOS, 
with an incompressibility of K = 380 MeV; and a = -356 MeV, (3 = 303 MeV and 7 = 
7/6 which corresponds to the so-called soft EOS with an incompressibility of K — 200 
MeV; C sym is the symmetry energy strength due to the difference of neutron and proton 
[ 26J, here C sym = 32 MeV is used to consider isospin effects, or C sym = for no isospin 
effect. V c is the Coulomb potential, JJ Yuk is Yukawa (surface) potential. In this model, the 
isospin effects can be included in some terms, such as in-medium nucleon-nucleon cross 
section and Pauli blocking [EHEEj- 

The time evolution of the colliding system is given by the generalized variational princi- 
pal. Since the QMD can naturally describe the fluctuation and correlation, we can study 
the nuclear clusters in the model [ EEl EZl CHI EHj In QMD model, nuclear clusters are 
usually recognized by a simple coalescence model: i.e. nucleons are considered to be part 
of a cluster if in the end at least one other nucleon is closer than Ar min < 3.5 fm in 
coordinate space and Ap min < 300 MeV/c in momentum space [ 26j. This mechanism 
has been extensively applied in transport theory for the cluster formation. 

3. Number-of-Nucleon Scaling of the Directed and Elliptic Flows 

Reactions of 40 Ca + 40 Ca, 86 Kr + 58 Ni and 86 Kr + 124 Sn at 25 MeV/nucleon have 
been simulated by IDQMD. To compare the results from different reactions, the similar 
centrality has been chosen for the above systems, namely 4-8 fm, 6-10 fm and 7-11 fm, 
respectively. Since the reaction systems tend to freeze-out around 120 fm/c in the model 
calculation, the physics results can be extracted in the freeze-out stage. In this work, the 
results are extracted at 200 fm/c. Some parts of physics results have been reported in 
our recent papers [ 30\ EH] • 

The directed flow Vi as a function of rapidity (y) for the above systems has been studied. 
Fig. [I] shows Vi versus rapidity (upper panels) and the nucleon- number scaled Vi/A as a 
function of rapidity. From this figure, we found that the slope of the directed flow as a 
function of rapidity is negative for three systems at 25 MeV/nucleon, which indicates that 
attractive mean filed is important in such a low energy [Elj. Before the nucleon- number 
scaling, the values of directed flow are different for different-mass clusters: the heavier 
the clusters, the large the absolute value of directed flow. However, all the curves of 
rapidity-dependent V\ almost collapse onto the same curves by dividing V\ by its number 
of nucleon, which illustrates that the directed flow of the light nuclear clusters satisfies 
the number-of-nucleon scaling. Actually, the previous experimental data for light nuclear 
clusters up to A = 4 for 86 Kr + 197 Au showed that the directed flow is approximately 
proportional to mass number [ . 

The upper panels of Fig. El show transverse momentum dependence of elliptic flow for 
mid-rapidity light fragments in three calculation cases: hard EOS without isospin related 
terms (hard_niso), soft EOS with isospin related terms (soft_iso), and soft EOS without 
isospin related terms (soft_niso). From the figure, it shows that elliptic flows are positive 
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Figure 1. Ca+Ca at 4-8 fm (left panels), Kr + Ni at 6-10 fm (middle panels), Kr + Sn at 7 - 10 
fm (right panels). The upper panels show v\ versus rapidity and the lower panels show wi/A, i.e. the 
number-of-nucleon scaled directed flow, versus rapidity. The dashed line guides the eyes. 



and increase with the increasing of p t even though the values are a little different for the 
different cases. The positive value of v 2 reflects that the light clusters are preferentially 
emitted within the reaction plane, and particles with higher transverse momentum tend 
to be strongly emitted within in-plane, i.e. stronger positive elliptic flow. In comparison 
to the elliptic flow at RHIC energies, the apparent behavior of elliptic flow versus p t 
looks similar, but the mechanism is obviously different. In intermediate energy domain, 
collective rotation is one of the main mechanisms to induce the positive elliptic flow [ 
EH 133 ESI E31 133 ESj In this case, the elliptic flow is mainly driven by the attractive mean 
field. However, the strong pressure which is built in early initial geometrical almond- type 
anisotropy due to the overlap zone between both colliding nuclei in coordinate space will 
rapidly transforms into the azimuthal anisotropy in momentum space at RHIC energies [ 
fT^j . In other words, the elliptic flow is mainly driven by the stronger outward pressure. 
The lower panels in Fig. [21 displays the elliptic flow per nucleon as a function of transverse 
momentum per nucleon, and it looks that there exists the number-of-nucleon scaling when 
Pt/A < 0.25 GeV/c. This behavior is apparently similar to the number-of-constituent- 
quarks scaling of elliptic flow versus transverse momentum per constituent quark (pt/n) 
for mesons and baryons which was observed at RHIC [Ej- 

4. ^4/^2 -scaling 

Recent the RHIC results show that v±jv\ for hadrons keeps an almost constant which 
is independent on p t [ 37j, we would like to know what the higher order momentum 
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Figure 2. The upper panels: Elliptic flow as a function of transverse momentum (p t ) for 86 Kr + 124 Sn 
at 25 MeV/nucleon. Three simulations were done: hard EOS without isospin related terms (hard_niso), 
soft EOS with isospin related terms (soft_iso), and soft EOS without isospin related terms (soft_niso). 
Squares represent for neutrons, circles for protons, triangles for fragments of A = 2, diamonds for A = 3 
and stars for A = 4; The lower panels: Elliptic flow per nucleon as a function of transverse momentum 
per nucleon. 
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Figure 3. The upper panels: v^/A as a function of pt/A for different particles, namely, neutrons 
(squares), protons (circles), fragments of A = 2 (triangles), A = 3 (diamonds) and A = 4 (stars) in three 
different calculation conditions. The middle panels: V4/A 2 as a function of (pt/A) 2 . The lower panels: 
the ratios of v^/v 2 for different particles vs p t . Three simulations were done in this figure as Fig. 2. 
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anisotropy in the intermediate energy is. In the present work, we explore the behavior 
of in the model calculation. Fig. [3| shows the features of V4 for 86 Kr + 124 Sn at 25 
MeV/nucleon for IDQMD simulation with hard_niso, soft_iso and soft_niso. Similar to 
the relationship of v 2 /A versus Pt/A ) we plot v^/A as a function of Pt/A (the upper panel). 
The divergence of the different curves between different particles in V4 j A versus Ptj A indi- 
cates no simple scaling of nucleon number for 4-th momentum anisotropy. However, if we 
plot V4/A 2 versus (p t /A) 2 (the middle panels), it looks that the points of different particles 
nearly merge together and it means a certain of scaling holds between two variables. Due 
to a nearly constant value of V4/V2 in the studied p t range (see the bottom panels in Fig. EI) 
together with the number- of- nucleon scaling behavior of v 2 /A vs ptjA, V4/A 2 should scale 
with (p t /A) 2 , as shown in the middle panels of Fig. El Interestingly, the constant value of 
v±jv\ does not depend on the EOS and isospin dependent term, which may indicate this 
constant value is robust for intermediate energy heavy ion collision. 

The RHIC experimental data demonstrated a scaling relationship between 2-nd flow 
(^2) and n-th flow (^ n ), namely v n (p t ) ~ v^ipt) [ 38j. It has been shown [EHIED! that such 
scaling relation follows from a naive quark coalescence model [E] that only allows quarks 
with equal momentum to form a hadron. Denoting the meson anisotropic flows by v n ^M(Pt) 
and baryon anisotropic flows by v n ^B{pt)) it was found that v^M(Pt) = (1/4)^2,m(p0 f° r 
mesons and v^siPt) = (1/3)^1 b(p*) f° r baryons if quarks have no higher-order anisotropic 
flows. Since mesons dominate the yield of charged particles in RHIC data, the smaller 
scaling factor of 1 /4 than the empirical value of about 1 indicates that higher-order quark 
anisotropic flows cannot be neglected. Including the latter contribution, one can show 
that 

^4,M 1 , lV4,q 1 1 V^ q 

7 r>~i _l _±_ 7 rv_/ __| _±_ 

v\ M ~ 4 2 v\ q ' v\ B ~ 3 3 v\ q ' 

where v n ^ q denotes the quark anisotropic flows. The meson anisotropic flows thus satisfy 
the scaling relations if the quark anisotropic flows also satisfy such relations. One can go 
one step further and assume that the observed scaling of the hadronic v 2 actually results 
from a similar scaling occurring at the partonic level. In this case, if one assumes [EH1EII1 
that the scaling relation for the partons is as follows: 

vl = K) 2 , (8) 

and then hadronic ratio v^/v^ then equals 1/4 + 1/2 = 3/4 for mesons and 1/3 + 1/3 = 2/3 
for baryons, respectively. 

If we assume the scaling laws of mesons (NCQ=2) and baryons (NCQ=3) (Eq. |7j) are 
also valid for A = 2 and 3 nuclear clusters, respectively, then V4/V2 for A = 2 and 3 clusters 
indeed give the same value of 1/2 as nucleons, as shown in Fig. Efc). Coincidentally the 
predicted value of the ratio of v^/v 2 for hadrons is also 1/2 if the matter produced in 
ultra-relativistic heavy ion collisions reaches to thermal equilibrium and its subsequent 
evolution follows the laws of ideal fluid dynamics [E]. It is interesting to note the same 
ratio was predicted in two different models at very different energies, which is of course 
worth to be further investigated in near future. 
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5. Summary 

In summary, we have investigated the behavior of anisotropic flows, namely v±, v 2 
and v 4 , for light nuclear fragments produced by 40 Ca + 40 Ca, 86 Kr + 58 Ni and 86 Kr + 
124 Sn at 25 MeV/nucleon for peripheral collisions in the framework of quantum molecular 
dynamics model. v 1 shows a negative slope versus rapidity, indicating the attractive 
mean field plays a dominant role in directed flow. Both i> 2 and V4 generally show positive 
values and increase with p t . The positive v 2 illustrates that the in-plane emission is 
preferential. By the scaling of the number-of-nucleons, both v±/A versus rapidity and 
v 2 /A versus Pt/A for light nuclear fragments approximately collapse on the similar curve, 
respectively, which means that there exists directed flow and elliptic flow scalings on the 
number-of-nucleon. Nucleon-number scaling originates from nucleonic coalescence. For 4- 
th momentum anisotropy i> 4 , it seems to be scaled by and v±jv\ is a constant, i.e. 1/2, 
for light nuclear fragments. Comparing our above predictions with intermediate energy 
HIC data for the flow scalings which is presented in this work is expected to shed light 
on nucleonic collectivity and its rich features. 
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